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Abstract. The extraction of the S-wave kaon-nucleon scattering lengths ao and ai from a combined analysis 
of existing kaonic hydrogen and synthetic deuterium data has been carried out within the framework of 
a low-energy effective field theory. It turns out that with the present DEAR central values for the kaonic 
hydrogen ground-state energy and width, a solution for ao and ai exists only in a restricted domain of input 
values for the kaon-deuteron scattering length. Consequently, measuring this scattering length imposes 
stringent constraints on the theoretical description of the kaon-deuteron interactions at low energies. 
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1 Introduction 

Recently, the DEAR collaboration at LNF-INFN has per- 
formed a measurement of the energy level shift and width 
of the kaonic hydrogen ground state [1] with a consider- 
ably better accuracy than the earlier KpX experiment at 
KEK [2]. The preliminary result of DEAR is 

e ls = 193 ± 37 (stat) ± 6 (syst) cV , 

r ls = 249 ± 111 (stat) ± 30 (syst) cV . (1) 

As can be seen from the above result, the accuracy is still 
tens of eV in the energy shift and more than 100 eV in the 
width. Now DEAR is being followed by the SIDDHARTA 
experiment that will feature new silicon drift detectors. 
The plans of the SIDDHARTA collaboration include the 
measurement of both the energy shift and the width of 
kaonic hydrogen with a precision of several eV, i.e. at the 
few percent level, by 2007. Moreover, SIDDHARTA will 
attempt the first ever measurement of the energy shift 
of the kaonic deuterium with a comparable accuracy and 
possibly, kaonic helium and sigmonic atoms. 

The necessity to perform measurements of the kaonic 
deuterium ground-state observables is justified by the fact 
that, unlike in the case of pionic atoms, the measurement 
of only the kaonic hydrogen spectrum does not allow - 
even in principle - to extract independently both S-wave 
KN scattering lengths ao and a\. This happens because 
the imaginary parts of these scattering lengths do not van- 
ish in the isospin limit, being determined by the decays 



into inelastic strong channels ttS, n°A, ■ ■ ■. Consequently, 
one attempts here to determine four independent quan- 
tities (real and imaginary parts of ao and ai) that re- 
quires performing four independent measurements - e.g., 
the energy level shifts and widths of kaonic hydrogen and 
kaonic deuterium. However, even though it is clear that 
ao and a\ can not be determined separately without mea- 
suring kaonic deuteron, it is still not evident whether it 
is possible to do so if one performs such a measurement. 
The reason is that the (complex) kaon-deuteron ampli- 
tude at threshold, which is directly determined from the 
experiment and which is expressed in terms of ao and a\ 
through the multiple-scattering series, is generally plagued 
by systematic uncertainties due to a poor knowledge of 
the low-energy kaon-nucleon dynamics. Thus, we have to 
understand in advance, whether these uncertainties arc 
small enough not to hinder a determination of ao and a\ 
from the forthcoming SIDDHARTA experiment. This is 
the main purpose of the present paper. 

The kaon-deuteron scattering process has been exten- 
sively studied in the past within potential (multiple scat- 
tering) approaches (see, e.g. [3-6]). To our knowledge, up 
to now the most detailed investigation of the problem 
within an effective field theory (EFT) framework is carried 
out in Ref. [7] on the basis of the so-called chiral unitary 
approach. In the present paper we undertake a new investi- 
gation of the problem within a systematic ficld-thcorctical 
approach which, at the lowest order, is formally similar to 
that of Ref. [7]. In addition to the previous studies, in this 
paper we also consider different sources of the systematic 
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error in the calculations, and try to give - whenever pos- 
sible - a crude estimate of the theoretical uncertainty. 

The main difference of the present article to previ- 
ous work is, however, that the existing approaches were 
exclusively concentrated on the prediction of the K~d 
scattering length from the input KN scattering lengths. 
We are not aware of the "reversed" analysis in the liter- 
ature, where the KN scattering lengths are determined 
from the input data of kaonic hydrogen and deuterium 
ground-state shift and width. However, this is exactly the 
type of the analysis that will be required in the near fu- 
ture for the SIDDHARTA data. In this paper - in the 
absence of any experimental data - we use "synthetic" 
data and show that the reversed calculations, owing to 
the non-linear dependence of the kaon-deuteron amplitude 
on the KN scattering lengths, turn out to be much more 
sensitive to the theoretical input on the deuteron struc- 
ture and the kaon-deuteron interactions, than a straight- 
forward evaluation of the K~d scattering length through 
the multiple-scattering series. This fact could potentially 
render a combined analysis of the hydrogen and deuterium 
data a beautiful testing ground for different EFT descrip- 
tions of the low-energy kaon-deuteron interactions and, 
as a result, might enable one to accurately determine the 
values of the scattering lengths a and a%. 

The paper is organized as follows. In section 2 we con- 
sider the extraction of the complex kaon-deuteron scat- 
tering length from the data on kaonic deuterium. The 
multiple-scattering series for the kaon-deuteron scattering 
amplitude is studied in section 3 and the issue of isospin 
breaking is addressed in section 4. Further, in section 5 
the numerical results of the simultaneous analysis of the 
existing kaonic hydrogen and synthetic kaonic deuterium 
data arc presented and discussed. Section 6 contains our 
conclusions. 



2 Kaonic deuterium 

In the experiments on hadronic atoms one measures the 
energy levels and widths of this sort of bound states. At 
present, there exists a well established systematic proce- 
dure for extracting the values of the pertinent hadronic 
scattering amplitudes at threshold from these measure- 
ments, based on non-rclativistic effective Lagrangians (see 
e.g. [8-13]). In particular, the case of the kaonic hydro- 
gen has been addressed in Ref. [12]. The calculation of 
the ground state energy of kaonic deuteron is completely 
analogous to the derivation of the Deser-type formula in 
the case of pionic deuterium, which was carried out in 
Ref. [13] - we do not repeat this derivation here. How- 
ever, as first pointed out in Ref. [14], in the case of kaonic 
atoms one should in addition check whether the relatively 
large decay widths of these atoms leads to a conflict with 
the use of the Raylcigh-Schrodingcr perturbation expan- 
sion for the energy levels. In that article it is proposed to 
estimate the quantity rEf s , where r denotes the lifetime 
and Ef s the binding energy for the ground state. Should it 
turn out that this product is of order one, using Rayleigh- 
Schrodinger perturbation theory would be questionable. 



The lifetime of kaonic deuterium is determined by the 
imaginary part of the K~d scattering length. In the ab- 
sence of any reliable experimental information, we use the 
crude estimate Im A^d — 1 fai which is consistent with 
the analysis of Ref. [15]. The decay width and the ground- 
state binding energy at lowest order in the fine-structure 
constant a are given by 

rf a = - ~ 4a 3 [4 lmA Kd ~ 1.2 keV , 
r 

Ef s ~ - n r a 2 ~ 10.4 keV . (2) 

This yields rEf s ~ 8.6, which is still large enough to jus- 
tify using perturbation theory (here, \x r denotes the re- 
duced mass of the kaon-deuteron system) . 

The (complex) kaon-deuteron scattering length A Kd 
can be extracted from the future SIDDHARTA data by 
using the Deser-type formula at next-to-leading order in 
isospin breaking, which is the same as in Ref. [13] 




= -2a 3 f4 A Kd 



x {l-2an r A Kd (lna-1) + •••}, (3) 

where the ellipses stand for small contributions which can 
be neglected at the accuracy we are working. On the other 
hand, since rEf s is not very large, it is not excluded that 
the corrections at next-to-next-to-leading order in isospin 
breaking (which are expected to amount to a few percent) 
should also be taken into account, when the accuracy of 
SIDDHARTA data is close to the planned one. It is how- 
ever, well-known that such calculations can be performed 
in a straightforward manner within the non-relativistic 
EFT and hence, no uncontrollable systematic error finally 
emerges at this place. One may therefore safely assume 
that the quantity Axd is directly determined from the ex- 
periment. 

3 Multiple-scattering series for the 
kaon-deuteron amplitude 

The main purpose of the present paper is to investigate 
whether the measurement of the quantity A^d, along with 
the K~p elastic scattering amplitude at threshold that is 
separately determined in the experiment on kaonic hy- 
drogen, enables one to extract precise values of ao and 
Oi. The first step is to express Axd in terms of ao and 
a± (and possibly, other physical parameters characteriz- 
ing the low-energy KN interaction) within the multiple- 
scattering theory. This procedure can be appropriately for- 
mulated by using the language of effective non-relativistic 
Lagrangians (see, e.g. [13, 16, 17]) that, in addition, allows 
one to systematically calculate the corrections. However, 
a strong predictive power of the non-relativistic approach 
is in fact based on a very subtle balance of different mo- 
mentum scales involved in the problem [13,22]. Namely, 
in this approach the couplings of the effective Lagrangian 
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are easily expressed in terms of the threshold parameters 
of the elementary KN interactions and the perturbative 
expansion of the kaon-deuteron scattering amplitude cor- 
responds to the usual multiple-scattering series. In anal- 
ogy with the case of the pion-deuteron system one would, 
however, expect that the three-body low-energy constant 
(LEC), which describes contact interaction between the 
kaon and two nucleons, is strongly enhanced and, since 
the value of this LEC is unknown, this leads to a large 
systematic error in the predicted kaon-dcuteron scatter- 
ing length [13, 16, 17]. On the other hand, if one uses the 
EFT with non-perturbative pions based on the Weinberg 
counting rules (see [24] for a recent review), both the di- 
mensional analysis and resonance saturation [13], as well 
as the study of the scale dependence [22, 23] carried out in 
the pion-deuteron case, indicate on a rather small uncer- 
tainty due to the three-body interactions - at the cost of 
the fact that now the expansion of the pion-deuteron scat- 
tering length should be carried out in Chiral Perturbation 
Theory (ChPT) and not in terms of the threshold param- 
eters of the underlying ttN scattering amplitudes. The lat- 
ter property may lead to even more serious problems in the 
kaon-deuteron case, owing to the non-perturbative charac- 
ter of the low-energy KN interactions and the necessity of 
using the chiral unitary approach (see, e.g. [18-21]). The 
solution to the above dilemma lies in the observation of 
a certain hierarchy of various contributions in the theory 
with non-perturbative pions, which could be described by 
the so-called modified power counting [22]. The existence 
of such a hierarchy, which can be traced back to the sup- 
pression of the diagrams with the creation/annihilation of 
virtual pions, indicates that the structure of the theory 
closely resembles that of the non-relativistic theory up to 
the three-body contact terms, whose value has to be deter- 
mined from the matching of these two theories. Further, 
from the direct comparison of the expressions of the pion- 
deuteron scattering lengths calculated within these two 
approaches one may conclude that the lowest-order three- 
body LEC in the non-relativistic theory can be effectively 
omitted, if the deuteron wave function, calculated in the 
EFT with non-perturbative pions, is used to evaluate ma- 
trix elements in the non-relativistic EFT. Below, we shall 
use this simple prescription for estimating the value of the 
three-body LEC in the case of the kaon-deuteron system 
as well. 

The pion-deuteron and kaon-deuteron systems differ 
in one crucial aspect. As it is well known (see, e.g. [7] 
and references therein), the multiple-scattering series for 
the kaon-deuteron scattering does not converge and re- 
quires a (partial) re-summation. This can be done most 
simply by using the so-called Fixed Center Approxima- 
tion (FCA), in which the nucleons are considered to be 
infinitely heavy. The validity of this approximation has 
been studied both in the potential scattering theory (see 
e.g. [25]) and in the EFT approach [26]. The fact that FCA 
can be a reasonable approximation even for Mk / fn p ~ 0.5 
(see Ref. [7] and references therein) is related to the pecu- 
liar cancellations at second order, which are discussed in 
Rcfs. [25, 26] (here, Mk and m p denote the masses of the 



charged kaon and the proton, respectively). In the second- 
order calculations in the non-relativistic EFT, which we 
have carried out, it is indeed possible to analytically verify 
the cancellation of the corrections to the FCA at leading 
order in Mfc/mp (here, one can use the corresponding 
analytic expressions from Ref. [13] after the substitution 
M n — ► Mk)- Moreover, it should be pointed out that in 
this theory the "binding corrections" [25] should be also 
included in the second-order term and, MS Mi result, the 
cancellation occurs in both isospin channels and not just 
in one as in Ref. [26]. Numerically, utilizing dimensional 
regularization, the corrections turn out to be of order of 
20 — 30% for the kaon-deuteron system and of order of 
a few percent in the pion-deuteron system. We expect 
these numbers to further decrease if the wave functions 
calculated in EFT with non-perturbative pions are used 
in the calculations. Note that, to the best of our knowl- 
edge, there exists no proof of cancellations beyond second 
order, but from e.g. the comparison to the exact solution 
of Faddeev equations [5] (see also the discussion in Ref. [7]) 
one may conclude that the numbers quoted above give a 
realistic estimate of the theoretical error due to FCA at 
all orders. It should however, be stressed that the EFT 
approach which is used in the present paper provides an 
appropriate tool for a systematic calculation of the correc- 
tions to the FCA, as well as the inclusion of higher-order 
(derivative) interactions. In order to provide a sufficient 
theoretical accuracy for analyzing SIDDHARTA results, 
such calculations might be necessary in the future. 

Using the FCA in the non-relativistic EFT and ne- 
glecting derivative interactions, we arrive at the expression 
for the kaon-deuteron scattering length, which is formally 
similar to the one from Ref. [7]: 



1 + th-f* 



(u 2 (r)+w 2 (r))a kd (r), (4) 



where M d is the deuteron mass, u(r) and w(r) denote 
the usual S— and D— wave components of the deuteron 
wave function, which are normalized via the condition 
J °° dr (u 2 (r) +w 2 (r)) = 1, and 



„ , . a p + a n + {2a p a n - b 2 x )/r - 2b 2 x a n /r 2 
akd(r) = ; . - , q h oa kd 



1 - a p a n /r 2 + b 2 x a n /r z 



with b 2 x = a 2 /(l + a u /r). Further, 



. M K 
1 H I a,- 



(5) 



(6) 



where a p ^ n ,x,u denote the threshold scattering amplitu- 
des for K~p — > K~p, K~n — > K~n, K~p — > K n and 
K n — > K n, respectively. Finally, the quantity 5a k d is 
proportional to the three-body LEC /q , which is a coun- 
terpart of the quantity /□ introduced in Ref. [13]. As dis- 
cussed above, one may assume that vanishes, if the 
deuteron wave functions u(r) and w(r) are those calcu- 
lated in the EFT with non-perturbative pions. We have 
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further carried out the dimensional estimate of the re- 
sulting uncertainty, which yields a few percent theoretical 
error in the kaon-deuteron scattering length. Thus, the 
whole procedure is consistent. 

The calculations with the deuteron wave functions of 
the non-relativistic EFT, which at the leading order are 
given by 

u(r) = y/ZyeTV, w(r) = , 7 2 ~ m p E d , (7) 

are more subtle (here Ed denotes the binding energy of 
the deuteron). Namely, expanding Eq. (5) in the multiple- 
scattering series corresponding to an expansion in pow- 
ers of 1/r and integrating this term-by-term, it is im- 
mediately seen that the subsequent terms diverge worse 
and worse as r — > 0. This means that the pertinent in- 
tegrals require renormalization and are scale-dependent 
(see e.g. Refs. [13,16,17]). On the other hand, the inte- 
gral before such an expansion is well-defined and does not 
have any scale dependence. It is however, clear that the 
uncertainty which is related to the existence of (large) 
three-body contact interactions can not simply disappear 
as a result of the re-summation of the multiple-scattering 
series. To solve this puzzle, note that the non-relativistic 
theory makes sense only for the momenta |p| <C A, with 
A being of the order of the pion mass. Hence, in order to 
effectively mimic the effect of the short-range physics, one 
could straightforwardly perform cutoff regularization of 
the multiple-scattering series and study the cutoff depen- 
dence of the calculated kaon-deuteron scattering length 
after the re-summation of the series (even now the limit 
A — > oo exists). This can be done most easily through the 
replacement 1/r — > {l — cxp(— Ar)}/r in Eq. (5). Fur- 
ther, in order to relate the value of A to the scale /z of 
dimensional regularization, which is usually used in the 
non-relativistic EFT, we calculate the mean value of the 
operator 1/r between the wave functions given by Eq. (7) 
twice: once using dimensional regularization with minimal 
subtraction and once using cutoff regularization. The com- 
parison of these two results gives A = sfe.[i ~ 1.65 \x (where 
e denotes the basis of natural logarithm) , so that the inter- 
val 100 MeV < fi < 250 MeV, which was fixed in Ref. [13], 
is mapped onto the interval 165 MeV < A < 412 MeV. 
Note also that, if instead of the wave functions Eq. (7) the 
wave functions from the EFT with non-perturbativc pions 
are used, the calculated kaon-deuteron scattering length 
has practically no A-dependence for 165 MeV < A < oo. 
This shows that it is consistent to use an "unregularized" 
1/r operator together with these wave functions, because 
the necessary cutoff is provided by the wave function itself. 



4 Isospin breaking 

The equation (5) contains four different combinations of 
the threshold amplitudes. Consequently, one has first to 
relate these amplitudes to the two scattering lengths ao 
and oi, which should then be determined from the anal- 
ysis of the combined data on kaonic hydrogen and deu- 
terium. In this work we take into account the leading-order 



Table 1. KN scattering lengths ao and ai (in fm) from the 
literature. These scattering lengths are used as an input in the 
calculations of the kaon-deuteron scattering length. 



Ref. 


a 


ai 


Meifiner and Oiler [19] 


-1.31+il.24 


0.26 + i0.66 


Borasoy et al, fit u [20] 


-1.48 + i0.86 


0.57 + i0.83 


Oiler et al, fit A4 [21] 


-1.23 + i0.45 


0.98 + i0.35 


Martin [27] 


-1.70 + i0.68 


0.37 + i0.60 



isospin-breaking corrections in the kaon-nucleon scatter- 
ing amplitudes which are due to the unitary cusps [12]. 
The re-summation of the bubble diagrams leads to the 
following simple parameterization 

_ I («o + ai) + Qoaoai _ 

a P ~ i I qo~l I -i ' a n — CLl, 

l + if(ao + ai) 

\ (op - ai) _ \ (q + ai) - ig e a ai 

QX ~ 1 (ao + ai) ' ' l-Hk(a + ai ) ' 

(8) 

where 

q c = s/2n c A , q = \Z2fi A , 

A = m n + Mjf o — m p — M K , 

m p M K m n M R0 , n . 

Mc = , , r , Mo = . , , (9) 

m p + M K m n + Mko 

with m n , being the masses of the neutron and the 

K°, respectively. Further, in the calculations we have used 
the input scattering lengths evaluated within various ver- 
sions of the so-called chiral unitary approach [19-21], as 
well as the experimental values from Ref. [27]. Table 1 col- 
lects these input values. Note that for us it is not possi- 
ble to straightforwardly use the scattering lengths (isospin 
basis) from Refs. [7, 18], since the definition of these scat- 
tering lengths (in the isospin limit) differs from the one, 
which is used here or in Refs. [19-21]. In particular, in 
Refs. [7,18] ao,ai are determined from the KN scatter- 
ing amplitude with CM energy taken at K~p threshold, 
whereas the masses of all particles in the isospin limit are 
taken equal to the average masses in the isospin multi- 
plets 1 . Note that, in contrary to this, the usual definition 
of the scattering lengths in the isospin limit of QCD plus 
QED implies that the CM energy is set exactly equal to 
the threshold energy. 

5 Numerical results and discussion 

The tables 2, 3 and 4 contain our results of the calculations 
of the kaon-deuteron scattering length with the use of the 

We are thankful to E. Oset and J. Oiler for a clarifying 
discussion on this and other related topics. 
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above formulae. In particular, in table 2 the results ob- 
tained by using the leading-order deuteron wave function 
in the non-relativistic EFT (see Eq. (7)) are displayed for 
different values of the cutoff parameter A (or, equivalcntly, 
//). Table 3 contains the results obtained by using the Paris 
and Bonn potential model wave functions as well as the 
deuteron wave functions calculated at NLO in the EFT 
with non-perturbative pions for two different values of the 
cutoff parameter A, which is present in this theory. From 
table 2 one sees - as expected from Rcfs. [f 3, 16, 17] - a 
rather strong dependence of the calculated kaon-deuteron 
scattering length on the cutoff parameter A, which has 
not disappeared after carrying out the re-summation of 
the multiple-scattering series. On the other hand, the de- 
pendence on the parameter A in table 3 is very weak (cf. 
with [22,23]) and all wave functions in this table yield 
practically the same result, which in most cases - except 
for the input from Ref. [21] - coincides with the result 
of table 2 somewhere around fi = M w . Consequently, this 
yields an useful "rule of the thumb" for a rough estimate of 
the value of the three-body LEC f$ in the non-relativistic 
EFT. 

Finally, in table 4 we show the effect of the isospin 
breaking on the theoretical value of the kaon-deuteron 
scattering length. Bearing in mind the large isospin-brea- 
king corrections in the KN amplitudes [12] as well as in 
the pionic deuterium [28], it is a little surprising that in 
most cases - again except the input from Ref. [21] - the 
leading-order isospin-brcaking effect at the end turns out 
to be very small. Note also that the large isospin-brcaking 
corrections, which are quoted in Ref. [7], most probably 
result from the peculiar definition of the scattering lengths 
in the isospin limit, which is adopted there. 

Next, we turn to the main goal of the present paper: 
how does one extract precise data of the kaon-nucleon 
scattering lengths ao and a\ from the combined measure- 
ments of the (complex) ground-state energy in the kaonic 
hydrogen and kaonic deuterium, which are carried out by 
the DEAR/SIDDHARTA collaboration? Let us start from 
the kaonic hydrogen case, for which the preliminary results 
of the measurements are already known. These measure- 
ments alone do not suffice to determine the values of both 
ao and oi, but yield a relation between these two quan- 
tities that leads to a strong restriction on the possible 
values that ao and ai can take. In order to understand 
how this restriction emerges, let us consider the expres- 
sion for the ground-state energy of the kaonic hydrogen at 
next-to-leading order in isospin breaking [12] 

eis - 7j As = -2«Vc a p {l - 2a/i c (lna- l)a p } .(10) 

With the help of this equation one may directly extract 
the quantity a p from the DEAR data. Using now Eq. (8) 
to relate the quantities a p and a , ai we obtain 

2go 2 «p n ni x 

a + ai + a ai - = . (11) 

1 - qoa p 1 - qoa p 

Together with the requirement Im aj > 0, which stems 
from unitarity, Eq. (11) defines a circle in the (Re a/, 




Re a z [fm] 



Fig. 1. Restrictions set by the DEAR data on the values of the 
scattering lengths ao and oi. For comparison, we give the scat- 
tering length calculations from different analyses: 1) MeiBncr 
and Oiler [19], 2) Borasoy, Nifiler and Weise, fit u [20], 3) Oiler, 
Prades and Verbeni, fit A4 [21], 4) Martin [27]. 



Ima/)-planc. Part of this circle is shown in Fig. 1 (note 
that, bearing in mind the preliminary character of the 
DEAR data [1], we use only central values in order to il- 
lustrate the construction of the plot and do not provide 
a full error analysis). In order to be consistent with the 
DEAR data, both ao and ai should be on the right of this 
universal DEAR circle. For comparison, on the same fig- 
ure we plot ao and a\ from table 1. As we see, in most of 
the approaches it is rather problematic to get a value for 
ao which is compatible with DEAR. This kind of analysis 
may prove useful in the near future, when the accuracy of 
the DEAR is increased that might stir efforts on the theo- 
retical side, aimed at a systematic quantitative description 
of the KN interactions within the unitarized ChPT. 

Further, from Eq. (11) one may determine e.g. a\. 
Substituting this expression into Eqs. (4), (5) and (8), 
one arrives at a non-linear equation for determining ao 
with a given input value of Axd- In the absence of exper- 
imental data on kaonic deuterium, we adopt the follow- 
ing strategy to solve this equation. In the literature, we 
find different theoretical calculations of the Kd scatter- 
ing length. We consider some of them, namely the results 
of Refs. [4,5,29,30], as synthetic experimental data and 
solve the corresponding equation with respect to ao- For 
completeness, we also present the results which are ob- 
tained from the same synthetic deuteron data together 
with the old KpX data on kaonic hydrogen, ei s = 323 cV, 
Ji s = 407 eV [2]. Note also that our solution includes only 
central values of the input kaon-deuteron and kaon-proton 
amplitudes. At this preliminary stage it is natural to post- 
pone the full error analysis, until the SfDDHARTA data 
on the kaonic deuterium become available. 

We start with the deuteron wave function in the non- 
relativistic EFT, which is given by Eq. (7). The results 
for the extracted scattering lengths ao and a\ arc given 
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Table 2. Kaon-deuteron scattering length And (in fm), calculated by using Eqs. (4,5) (setting Sa,Kd = 0) and the input from 
table 1. The deuteron is described by the wave function in Eq. (7). The values of the regularization parameter A shown in this 
table, correspond to the following values of the dimensional regularization scale: fi = 100 MeV; /j, — \a — 250 MeV; [i = oo 
(cf. with Ref. [13]). 



Ref. 


A = 165 MeV 


A = 231 MeV 


A = 412 MeV 


A = oo 


[19] 


-1.35 + i 1.58 


-1.31 +i 1.42 


-1.21 + i 1.24 


-1.07 + i 1.13 


[20] 


-1.63 + i 1.39 


-1.56 + i 1.24 


-1.41 +i 1.08 


-1.25 + i 0.99 


[21] 


-1.30 + il.08 


-1.30 +i 0.96 


-1.18 + i0.82 


-1.01 + i 0.75 


[27] 


-1.70 + i 1.26 


-1.61 +i 1.14 


-1.46 + i0.99 


-1.31 +i 0.90 



Table 3. The same as in table 2, but with the deuteron wave functions calculated in the Paris and Bonn potential models 
as well as at NLO in EFT with non-perturbative pions for two different values of the cutoff parameter: NLO(l) stands for 
A = 450 MeV and NLO(2) for A = 600 MeV. The parameter A is set to oo. 



Ref. 


Paris 


Bonn 


NLO (1) 


NLO (2) 


[19] 


-1.29 +i 1.66 


-1.28+il.65 


-1.29 + i 1.66 


-1.30 + 11.67 


[20] 


-1.55 + il.48 


-1.54 + i 1.47 


-1.55 + i 1.48 


-1.56 + il.49 


[21] 


-1.10 + il.l0 


-1.10 + i 1.09 


-1.09 + i 1.10 


-1.11 + i 1.11 


[27] 


-1.66 + i 1.28 


-1.65 + i 1.27 


-1.66 + i 1.28 


-1.67 + i 1.29 



Table 4. Isospin-breaking corrections to the Kd scattering length for different input values of ao and a\. NLO EFT wave 
functions are used for the deuteron. The label "sym" refers to the scattering length evaluated in the isospin symmetry limit. 



Ref. 


NLO(l), sym 


NLO (1) 


NLO(2), sym 


NLO(2) 


[19] 


-1.26 +i 1.65 


-1.29 +i 1.66 


-1.27 + i 1.66 


-1.30 + 11.67 


[20] 


-1.48 + i 1.46 


-1.55 + i 1.48 


-1.50 + i 1.47 


-1.56 + i 1.49 


[21] 


-0.85 + i 1.00 


-1.09 + i 1.10 


-0.87 + i 1.01 


— 1.11 + i 1.11 


[27] 


-1.64 + i 1.22 


-1.66 + i 1.28 


-1.66 + i 1.24 


-1.67 + i 1.29 



in table 5, where we had to choose a rather small value 
of the cutoff parameter A = 132 MeV corresponding to 
\i = 80 MeV. From this table one already sees the main 
property of the solutions: due to a highly non-linear form 
of the equation which determines ao and a\ , the solutions 
do not always exist (at least for those values of these scat- 
tering lengths which can be termed physically reasonable) . 
Note that this mainly concerns the solutions which arc 
obtained using DEAR input. The KpX data turn out to 
be much less restrictive. For example, if we increase the 
value of the parameter /j, even up to 100 MeV, the sole 
solution with DEAR input in table 5 disappears, whereas 
the solutions with KpX input are still present. Note also 
that we have analyzed more synthetic input than is fi- 
nally shown in table 5: for most of the input values of the 
kaon-deuteron scattering length which are known in the 
literature, there exists no solution for ao and a±. 

Bearing in mind the results in table 5, it is not sur- 
prising that, using the wave functions calculated an NLO 



in the theory with non-perturbative pions, the solutions 
exist only with the KpX input. The corresponding results 
are shown in table 6. 

One might finally ask the following question: how large 
should the kaon-deuteron scattering length be so that a 
solution for ao and ai exists at all? In order to answer 
this question, we have scanned the (ReAxd, Im Air- 
plane in the interval —2 fm < KeAxd < and 0.5 fm < 
Im And < 2.5 fm and tried to find solutions, using DEAR 
input data. The results of this investigation, which are 
displayed in Fig. 2, are very interesting: it turns out that 
the solutions exist only if ImAxd 1 fin and moreover, 
if Im Axd — 1 fm then one finds solutions only in a very 
small interval around KeAxd — —1 fm. Some representa- 
tive solutions are shown in table 7. We wish to also note 
that all this agrees with the scattering data analysis, car- 
ried out in Ref. [15]. If Imixd crosses the border of the 
shaded area in Fig. 2 continuously from below, then on the 
same branch one gets the solution with Imai < that is 
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Table 5. The scattering lengths ao and ax, obtained from a simultaneous analysis of DEAR (KpX) data and the synthetic input 
for the Kd scattering length. The deuteron wave function is described by Eq. (7) and the regularization parameter A = 132 MeV. 



Input A K d 


DEAR 


KpX 


-1.34 + i 1.04 [4] 


ao = -1.27 + i 0.36 
oi = 1.19 + i0.09 


a = -1.34 + i 0.63 
ai = 0.01 +i0.55 


-0.85 + i 1.10 [5], Faddeev Eq. 




a = -1.40 + i 0.82 
ai = -0.09 +i 0.07 


-0.75 + i 1.12 [5], FCA 




a = -1.41 +i0.84 
ai = -0.10 + i 0.01 


-0.78 + i 1.23 [29] 




a = -1.43 + i 0.82 
ai = — 0.04 + i 0.01 


-1.92 + i 1.58 [30] 




a = -1.64 + i 0.63 
ai = 1.03 + i 0.02 



Table 6. The same as in table 5, but with the NLO EFT wave functions and A = oo. Only KpX input is displayed, since there 
are no solutions with DEAR input. 



Input Axd 


A = 450 MeV 


A = 600 MeV 


-1.34 + i 1.04 [4] 


a = -1.36 +i 0.59 
a x = 0.10 + i 0.63 


a = -1.35 + i 0.60 
ax = 0.07 + i 0.62 


-0.85 + i 1.10 [5], Faddeev Eq. 


a = -1.41+i0.80 
a x = -0.04 + i 0.08 


a = -1.41 + i 0.80 
oi = -0.05 +i 0.08 


-0.75 + i 1.12 [5], FCA 


a = -1.42 + i0.82 
ax = -0.04 + i 0.03 


a = -1.42 + i 0.82 
oi = -0.05 +i 0.03 


-0.78 + i 1.23 [29] 


ao = -1.44 + i 0.80 
ax = 0.01 + i 0.03 


ao = -1.44 + i 0.80 
ax = 0.00 + i 0.03 



forbidden by unitarity. On the other hand, if KpX input 
data are used, the solutions exist in a much larger area. 



6 Conclusions 

The message of our investigation is very clear. Up to now, 
in the theoretical description of the deuteron we have re- 
stricted ourselves to the lowest-order expression in the 
non-relativistic EFT (derivative interactions are neglec- 
ted) and, in addition, worked in the limit of infinitely 
heavy nucleons (FCA). It is widely believed that these 
approximations represent a good starting point for the 
description of the kaon-deuteron scattering length. More- 
over, in the non-relativistic EFT the corrections are sys- 
tematically calculable. 

Within the approximations described above it turns 
out that the combined analysis of DEAR/SIDDHARTA 
data on kaonic hydrogen and deuterium is more restric- 



tive than one would a priori expect. In particular, we 
see that solutions exist only in a rather small area of the 
(Re Axd, Im^4if £ ;)-plane. Due to this fact, in certain cases 
it might be possible to pin down the values of ao and a\ at 
a reasonable accuracy, even if Axd itself is not measured 
very accurately. Moreover, if the corrections to the lowest- 
order approximate result, as expected, are moderate, they 
should not change the qualitative picture. In our opinion, 
it could prove very useful to perform the calculations of 
these corrections before starting to analyze the forthcom- 
ing SIDDHARTA data since, as one may conclude from 
the discussion in the present paper, exactly these correc- 
tions constitute the largest potential source of theoreti- 
cal uncertainty at present. On the other hand, crude esti- 
mates show that the uncertainty, coming from the short- 
range QCD dynamics (three-body LEC) is rather small 
and should not likely hinder the analysis of the future SID- 
DHARTA data on the kaonic deuteron. Note also that the 
above conclusion concerns DEAR input only. The analysis 
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Fig. 2. The region in the (ReA^d, ImAxd)~plane where so- 
lutions for ao and ai exists. The cutoff parameter is A = 
600 MeV. For comparison, we also show the results of cal- 
culations of the kaon-deuteron scattering length from table 3 
(filled diamonds), the synthetic data from tables 5, 6 (squares) 
and our representative solutions from table 7 (circles). 

Table 7. Representative solutions for ao and ai using the NLO 
EFT wave functions together with DEAR input and Re Ak<i — 
— 1 fm, ImAftTd ~ 1 fm. 



A K d 


A = 450 MeV 


A = 600 MeV 


-1.10 + i 1.00 


a = -1.25 + i 0.38 
ai = 1.06 +i 0.00 


a = -1.24 + i 0.38 
ai = 1.02 + i 0.03 


-0.80 + i 0.90 


a = -1.07 + i 0.43 
ai = 0.44 + i 0.07 


a = -1.06 + i 0.42 
oi = 0.42 + i 0.08 



using KpX input turns out to be much less restrictive and 
therefore less informative than with the DEAR input. 

To summarize, one may expect that the combined 
analysis of the forthcoming high-precision data from 
DEAR/SIDDHARTA collaboration on kaonic hydrogen 
and deuterium will enable one to perform a stringent 
test of the framework used to describe low-energy kaon- 
deuteron scattering, as well as to extract the values of 
ao and a\ with a reasonable accuracy. However, in or- 
der to do so, much theoretical work related to the sys- 
tematic calculation of higher-order corrections within the 
non-relativistic EFT is still to be carried out. 
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